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Friction velocity u*and roughness length z 0 of atmospheric surface 

boundary layer in sparse-tree land 1 
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Abstract Sparse-tree land is one of the typical lands and can be considered as one typical rough 
surface in boundary layer meteorology. Many lands can be classified into the kind surface in the view 
of scale and distribution feature of the roughness elements such as agroforest, scatter planted or 
growing trees, savanna and so on. The structure of surface boundary layer in sparse-tree land is analyzed 
and the parameters, friction velocity u* and roughness length z 0 are deduced based on energy balance law and 
other physical hypothesis. The models agree well with data of wind tunnel experiments and field measurements. 
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Introduction 

Protection forest is one of the most frequently used 
measures to protect crops from strange wind and to 
reduce soil evaporation and soil erosion. Protection 
forest in farmland can be classified in three kinds 
(Cao 1983): One kind is shelterbelts. Second is 
sparse tree in farmland, generally called agroforest in 
China. This kind includes regularly and randomly 
planted or naturally growing trees in farming area. 
The third is cluster of trees. These three kinds of 
forest exist together in farmland area. The shelter- 
belts are general and the protective effect and aero¬ 
dynamics were very deeply studied (Bradly 1983; 
Cao 1983; Eimern 1963; Jensen 1961; Plate 1971; 
Schwartz 1995; Wei 1987; Wang 1996; Zhang 1984; 
Zhu 1993). But the researches of aerodynamics of 
second and third were few. Only some micrometeo- 
rological measurements were reported (Li 1994; Lu 
1983). Some related studies aimed at the surfaces 
with solid rough elements (Mashall 1971; Raupach 
1980; Wooding 1973) and the physical characteris¬ 
tics of these elements are different from trees. It is 
important to know the relationship between the aero¬ 
dynamics of the surface and the structure of 
sparse-trees because this is not only necessary to 
understand protection effect of the forest, but also 
helpful to provide parameters to climatic models. 
This paper will deduce the models of friction velocity 
and roughness length of surface boundary layer in 
sparse-tree land and demonstrate the models with 
data of wind tunnel experiments and field measure¬ 
ments. 
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Structure of surface boundary layer in 
sparse-tree land 

In neutral stratification, surface wind velocity profile 
in open land is 

u 0 (z)=— In— (1) 

K Zoo 

Where u. 0 , * =0.4, z 00 are friction velocity, Karman 
constant and roughness length of the open surface 
respectively. 

When the surface flow approached the sparse-tree 
land, the surface boundary layer change its structure 
and become relative stable, as the distance is far 
enough. The boundary layer can be roughly divided 
in three layers as follow (Fig.1): 

First layer is the boundary layer immediately on 
the ground that is generated by the friction between 
airflow and the ground. This layer is very thin and its 
ceiling is under the tree height. The airflow of the 
layer is not uniform in horizontal because of the dis¬ 
turbance of the trees. But except for the area very 
close to trees, the wind velocity profile in the layer 
can be approximately expressed in logarithm 

«'(*)-— In— (2) 

K Z ' Q 

Where u* and z 0 ’ are friction velocity and roughness 
length of the surface respectively. And generally u*< 
u*o. We call this layer “under layer”. 

Second layer is on “under layer”, ranging from 
ceiling of “under layer” to the height about 2H. This 
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layer is generated mainly by the friction between the 
air flow and the trees. The thickness of the layer is 
depended on stratification. The parameters such as 
friction velocity and roughness length are also influ¬ 
enced by the structure of the trees (canopy density 
and individual tree’s porosity or permeability). We call 
this layer “roughness sublayer”. The wind profile in 
the layer is 


width of tree respectively; r is relative approaching 
wind velocity to the tree; u e 0 is geometrical mean 
wind velocity in open land integrated from z 0 to H. It 
was approximately expressed as (Zhu 1993) 


k ez m 


u(z) =—+<p) (3) 

k z 0 

Where u. and z Q are friction velocity and roughness 
length of the layer respectively. <£is influence term to 
indicate the difference between logarithm profile and 
actual profile, d is zero-plane displacement, which is 
function of canopy density C of the trees (Guan, 
1997) 


d = 


0 

H 

a 


In— 


C ^0.116 
C> 0.116 


(4) 


where a=2.735 is an empirical constant. C 0 =0.116 is 
the maximum C for d=0 

The third layer is on the “roughness sublayer”, The 
wind profile in the layer is 

. . u, z —d 

u(z) =— In- (5) 

K z 0 

Where all the symbols are same as equation (4). 
This layer is called “inertial sublayer” 

We call “roughness sub layer” and “inertial 
sublayer” together as “internal boundary layer” 
(IBL). 


Friction velocity u* in IBL 

According to the principle of constant flux in sur¬ 
face boundary layer, the momentum flux M down to 
the surface in IBL is equal to the sum of drag of the 
trees x ,and the stress of the ground x 2 which are 
analyzed follow. 

Drag of the individual tree is 

t x =^C T pHD(ru, 0 ) 2 (6) 


and C T is drag coefficient of isolated tree which is 
expressed as (Guan, 1999) 

C T =— A(l-a 2 ) (8) 

r 2 

Where A=0.92 is empirical constant and a is perme¬ 
ability of isolated tree. 

If the trees are considered as cylinder shape and 
tree density is one on ground area s, then 

s= h (D/2) 2 /C (9) 

considering (6) to (9), drag of trees on unit ground 
area is 


r i= r i'/s= 

p CAH(1 - a 2 ) (rUeo) 2 /( x D) (10) 

The stress on ground is approximately expressed as 

t 2 = p (ru*o) 2 (11) 

Based on the wind tunnel experiment and field 
measurement data the empirical relation of r and C 
was 

r =0.88e' 40C + 0.12 (12) 

with correlation 0.986. 

Considering 

M= r,+ r 2 (13) 

M = pu* 2 (14) 

and (6) to (11), then 


w. 


— = (0.88e'* c + 


M* n 


0 . 12 ) 


.2 \ 1 H 2 


, 1/2 


CAHQ. - a )——z (In- Y +1 


7IK 


ez 


oo 


(15) 


where p is air density; H and D are height and crown 
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Model values and measurement in wind tunnel ex- pared in Table 1. The model agrees well with the 
periments and field of u*/u*o (Guan, 1997) are com- measurements. 



Fig. 1. Sketch map of structure of boundary layer in spars e-tree land 
(vertical distance not drawn to scale) 


Table 1. Ratio of friction velocity u*/u«o in boundary layer above sparse trees to open land 


Canopy density 


Wind tunnel experiment 


Land measurement 

0.056 

0.113 

0.225 

0.450 

0.900 

0.60 

0.80 

Measurement values 

2.92 

2.96 

3.10 

2.33 

2.03 

2.06 

1.98 

Model values 

2.81 

3.15 

2.94 

2.09 

1.82 

1.88 

1.80 

Errors 

0.11 

-0.09 

0.16 

0.24 

0.21 

0.18 

0.18 


Roughness length of IBL 

From (3) we deduce 


with correlation 0.980. 
Let 


ia 

H 


(1 —~) ex P(^ - ) 

H u. 


a = f}\n{H/zoo)= 

( 16 ) (0.55e 7OC +0.45)ln(H/zoo) 


If the approaching wind velocity relative to open 
land at height H is n, then 

(it) 

K Zoo 

So that (16) becomes 


then 


Zq 

H 


= (1—~)exp(^-7— 

H u. /«.„ 


) 


( 20 ) 


From measurement data and calculation by model 
(15) we deduced that <£is related to C as follow 


is 

H 


(!"—) exp 
li 


u. // 

—(7'In—) 


u 


*0 


‘•00 


(18) 


J0.9 c ^0.2 

y45C C < 0.2 


( 21 ) 


in which all terms are known or can be evaluated 
except for 4and <t>. Based on wind tunnel data, n 
was depended on the canopy density of sparse trees 
and the empirical relation was 

n = 0.55e 70C + 0.45 (19) 


with correlation 0.966. 

The model values and the “measurements” are 
shown in Fig.2. The “measurements" of ZqIH were 
the results of calculation by approximation method 
based on wind tunnel and field data (Guan 1997). It 
can be seen that the model agrees well with “meas¬ 
urements”. 
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Fig. 2. Comparison of the model values (line) and the 
“measurements”(dots) 


Conclusions and discussions 

Sparse-tree land is one of the typical lands and 
can be considered as one typical rough surface in 
boundary layer meteorology. Many lands can be 
classified into the kind surface in the view of scale 
and distribution feature of the roughness elements 
such as agroforest, scatter planted or growing trees, 
savanna and so on. This paper gives the models of 
friction velocity u. and roughness length z 0 of this 
kind surface based on physical laws and some em¬ 
pirical relations. The models agree well with meas¬ 
urements and calculated results from measurements. 

The authors think that this study is preliminary and 
more deep analysis and measurement are expected. 
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